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] Proceeded from the gravitation equations proposed by one of authors it was argued in a previous paper that there can exist 
O |, supermassive compact configurations of degenerated Fermi-gas without events horizon. In the present paper we consider the 
' stability of these objects by method like the one used in the theory of stellar structure. It is shown that the configurations 
with an adiabatic equation of state with the power > 4/3 are stable. 
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J_, 1. Introduction 

Thirring [Thirring 1961 proposed that gravitation could be described as a tensor field tpafiix) of spin two in 4 



X 



dimensional Pseudo- Euclidean space-time i?4, where the Lagrangian, describing the motion of test particles in a 
^) , given field, is of the form 

L = -™pc[5a/3Wi"il'^' • (1) 

In this equation g^p is a tensor function of il^ap, "nip is the mass of the particle, c is the speed of light and 
^ , i" = dx" /dt (Greek indices run from to 3). 

A theory based on such a Lagrangian must be invariant under some gauge transformations ipap — 4'af3 that 
are a consequence of the existence of " extra" components of the tensor tpap ■ Transformations tpaf3 — > "00/3 give 
■ rise to transformations gap — > gap- Therefore, the field equations for gap{x) and equations of the motion of a test 
particle must be invariant under these transformations of the tensor g^^. 

The equations of motion of the test particle, resulting from eq. ([T]) are also equations of geodesic lines of 
the 4-dimensional Riemannian space-time Vn whose the metric tensor is gafsiip)- Therefore, if the transformation 
TpafS — ^ V'a/3 leaves the equation of the motion invariant, then the corresponding transformation g^p — *■ gaf3 is 
some mapping V — > V of the Riemannian spaces leaving geodesic lines invariant, i.e. a geodesic (projective) 
one. If not only eq. ([T]) but also the equations of field contain ipap in the form gapiip), then clear that only 
geodesic-invariant equations for gap are permissible in such a theory. Such of kind equations which have not 
physical singularity in the spherically-symmetric field was proposed in the paper [Verozub 1991] . These equations 
are of the form 

Kpn - ^'-sBi = 0. (2) 
The equations are vacuum bimetric equations for the tensor 



Kp^Kp-^lp^ (3) 



where 

Kp = Kp + 1)-' [sZK^ + sjri^] , (4) 
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are the ChristofFel symbols of space-time E4 whose fundamental tensor is -qap, F^,^ are the Christoffel symbols 
of the Riemannian space-time V4, whose fundamental tensor is ga/i and 6'^ is the Kroneker delta. The semi-colon 
in eq. ^ denotes the covariant differentiation in £'4 0. 

Eqs. ^ are invariant under arbitrary transformations of the tensor gaf3 retaining invariant the equations of 
motion of a test particle. Thus, the tensor field gaf3 is defined up to geodesic mappings of space-time V4 (in the 
analogous way as the potential Aa in electrodynamics is determined up to gauge transformations). A physical 
sense in the theory have only geodesic invariant values. The simplest object of that kind is the object B'^^ which 
can be named the strength tensor of gravitation field. The coordinate system is defined by the used measurement 
instruments and in each case is a given. 

Because of the gauge (i.e. geodesic or projective) invariance, additional conditions can be imposed on the tensor 
gafi- In particular [Verozub 1991] . under the condition 

Qc - - rL - (7) 

eqs. ^ are equivalent to the following system: 

Raf3 — 0, (8) 
= 0, (9) 

where Rap is the Ricci tensor. 

Proceeded from these equations in the papers [Verozub 1996] - [Verozub 1997] . it was shown that there can 
exist supermassive equilibrium compact objects without events horizon. Maybe just such of kind objects are in the 
galactic centers [Verozub 1999] . 

A geodesic-invariant generalization of the equations ([2|) inside matter can be found in the following way (see 
also [Verozub and Kochetov 2000] ). 

Transformations of the Christoffel symbols under the geodetic (i.e. projective) mappings is of the form: 

fj^ = r^^ + ^pS!^ + ^,6-^, (10) 

where tp-y is a vector- function of x". This equation has a simple interpretation in an 5- dimensional manifold M5, 
where the admissible coordinates transformations arc given by 

=i''{x",x\x^,x^), (11) 

X* =x^- J ific^dx". (12) 

Namely, eqs. (fTOj) can be interpreted as the transformation of 4-components F^^ of the connection coefficient Tg^ 
{A, B,C = 0..4) in Ai^ under the transformation (fT2|) . if the components F"^ obey the condition F"^ — (5^ . 
For this reason we will consider the tensor g^p as 4-components of n + 1- dimensional tensor 



54a 544 

The components gap are transformed under (|12p as follows: 

5a/3 = 9af3 + QiiV'a^lS + dAa^f) + Qifi^a, (14) 
5a4 = 9ai + gMfa, (15) 
544 = 544- (16) 

Transformations of the components F^^ under the geodesic mappings are given by 



^There was a misprint in eqs. (9) of the paper | Verozub 19961 . These equations must be read as follows: 

limA = l, hm{B/r^) = l, lim C = 1. (6) 
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where n is dimension of our space-time. Therefore, such a transformation for Qa coincides with (|15p . if 344 = n + 1. 
For this reason, we wiU assume that 

Then there exists the geodesic-invariant tensor 

Gap ^ gap ~ {n+iy^QaQp, (18) 

and the geodesic-invariant generahzation of the Einstein equations with matter source are of the form 

Kp;,- BZaB^Py = k {Tap - l/2GapT) , (19) 

where k — 8ttG/c^, G is the gravitational constant, c is speed of light, Tap is the matter energy-momentum tensor, 
T = G^^Tap. Thus, we assume that inside matter the gravitation equations under consideration just as the vacuum 
equations coincide with the Einstein equations at the gauge conditions Qa = 0. 

In the previous paper [Verozub 1996| the objects was considered as homogeneous. In the present paper we find 
the solution of the gravitation equations inside the objects and obtain more rigorous prove of the objects stability. 



2. The Internal Structure 

Consider here the relativistic equations of the objects structure. In the spherically symmetric field the metric 
differential form of space-time V4 is given by [Verozub 1991] 

ds^ = C dx°^ -Adr'^-B [dO"^ + sin^ dip"^] (20) 

where A, B and C are the functions of the radial coordinate r. 

Temporarily let us replace in eq. (|20p the radial coordinate r by / = \/B: 

ds^ = dx"^ - e'' df - f{de^ + sm^Od'^}^), (21) 

where 

df/dr 



= C, e'' = -4^. (22) 



We must solve the system of the equations 



For an ideal fluid eqs. ((23)) are 



Rap = k{Tap - -^gapT), (23) 
Qa = 0. (24) 



: kp, 

kfi, (25) 





' p) 


1 






^P 
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1 
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^? 






to 


2 


/ 



where Aj = dX/df, 7^ = dj/df, X'jj — d^\/dp, = pp, p is the matter density and p is the pressure. 
These equations are equivalent to the following system [Tolman 1969| : 



kp, 
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2 = kfi, (26) 



df 

The general solution of this system is given by 



e> = a,(l-^)e*, (27) 



l-^y\ (28) 



/ 

dp _ 1 (m + p) + + kpp 



(29) 



df 2/2 1 - (* + ^72)// ' 

where 

* = A; / (30) 



' l-(* + (72)//' ^ ^ 



(Ti and tT2 are constants. Returning to the variable r we obtain finally three equations: 

C = .r (l-^^e* 



c?p ^ 1 (/x+p) ('? + a2 + fcp/^) rf/ 
dr ~ 2 p 1- (* + a2)// dr' 

The gauge condition ([7]) yields 4th equation: 

ACA-4=a3, (33) 

where 173 is a constant. 

The eqs. (f32l) and ([33]) are the general solution of the used gravitation equations inside and outside objects 
under consideration. 

The constant because of the conditions ^ is equal to 1. The constants CTi and (J2 for the solution outside 
the considered objects can be found from the conditions at the infinity Setting in eq. I|32p ^ = and p = we 
obtain equations 

C = a,[l-jy (34) 

which together with ([33]) yields the solution of the vacuum gravitation equations ([2]) [Verozub 1991] . The constants 
are: ax = 1, a2 = rg = IGMjc? , where M is the mass of the considered object. 

Consider the solution of eqs. ([32]) inside the objects. From the condition that at r = i?, where R is the radius 
of the object, the right hands of eqs. ([34]) and (|35|) must coincide with the ones of the first two eqs. ([32]) we find 

CTi = exp[-^'|^^^(fl)], (36) 

(72 =rg-$|^^^(fl). (37) 

By using the notions 

u = - *, (38) 

= exn . ,,9^ 
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we obtain finally the following system of the equation which determine the internal structure of the object: 



dr f 



1 = (40) 



dp ^ 1 {^J,+ p)r'^ kpp + rg-u 

dr 2 P f-rg+u ^^'^ 

— = -kfir^v, (42) 
dr 



dv k {fi + p)r 



2 



v' (43) 



dr 2 / - rg 
and the boundary conditions: 

/(i?) = (i?-Vr3)i/3, p{R) = 0, uiR) = 0, viR)^l. (44) 



The function v very little differs from 1 because the power of the exponent (|39p is much less than 1. Therefore, 
it follows from the ((40|l that the function / very little differ from the ones in the solution of the vacuum equations 
@. 

At i? 3> the function / ~ r. At this condition these equations coincide with the ones of general relativity. 
We use an approximation equation by Harrison [Harrison 1965j for the baryons density as a function of the 
matter density p which is right from 8 to at least 10"'^^ g/cm? : 

n6 = gip(l + 92p9/i6)-4/9^ (45) 
where qi = 6.0228 • lO^^ and 92 = 7.7483 • 10"^" in CGS imits, and the equation of the state 

p) c'. (46) 



drib/ dp 



In addition to the ordinary solution of eq. (pUl) (i. e. configurations of the white drafts and neutron stars) there 
exist solutions with large masses from IO^Mq up to 10"'^'' Mq. Fig. [1] show the distribution of the density p{r) 
inside the configuration with the mass 2.6 • lO^M© (its radius is equal to O.O57i?0). 

Figs. [21 [3] and m show the relations " central density- mass" , "central density-radius" and "mass-radius" for 
the objects under consideration. 

There are no solutions for greater masses. 

In according to figure [4] the relation " mass-radius" for these configuration is given by : 



^ = 4.07.10-(-^)- . (47, 
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Fig. 2: The relation between the cen- 
tral dencity and the object radius 



Fig. 3: The relation between the cen- 
tral density and the object mass 



^ Fig. 4: The relation between the 

"O"^ masses and radiuses of the objects 
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3. Stability of Equilibrium Configurations 

It is not clear, at present, whether the found solutions for the big masses make up a continuous class of solutions with 
the solutions for neutron stars and dwarfs or not. Besides, the above figures have not any extremes points. For this 
reason the figures do not give us direct evidences of the stability or instability of the objects under consideration. 

The stability of the objects was argued first in paper [Verozub 1996] where the objects was considered as 
homogeneous. Following a classical method by [Ledoux 1958] and [Cox 1980] consider now the problem of the 
stabihty of the configurations more rigorously, without the assumption of homogeneity. 

The complete energy of a spherically symmetric object of the mass M and the radius R can be written as 
follows 

M 

E= (w + 0)(im, (48) 



where m is the mass of the matter inside the sphere of the radius r, u is the intrinsic energy, (j) is the gravitational 
potential. The function C ^ u + (f> will be considered as the function of the variables m, r(m) and r'(m) = dr/dm. 
Consider the functional 

M 

E = J £[m,r{m),r'{m)]dm (49) 



on the set of all continuously differentiable function r(m) with the boundary conditions r(0) = 0, r{M) — R and 
find at which conditions the first variation SE is equal to zero and the second one is a positive at small isoentropic 
disturbances r — > r + 6r. 
Since 

the condition 5E = yields the Lagrange equation and end condition. 
The Lagrange equation 

dm \dr' J dr 

is given by 

_ (7 r _ / r)-?; \ 

-.9 = 0, (52) 



d 

dm 



du 
dp 



si 



where {d/dp)g denotes the derivative at a constant value of the entropy S, g — F/m^ and the force F affecting 
the particle with mass mp is given by [Verozub 1996] 



where / = (fg + r^)^/^. 
Since the pressure is 



F = -^II^[l-'j], (53) 



dp 



the eq. is the equation of hydrostatic equilibrium 



dp 
dr 



pg. (55) 
When the Lagrange equation (|5T]) is satisfied, the end condition can be obtained from the equality 



M 

dC 



= 0, (56) 





where rj = Sr. Because of spherical symmetry we have rj (0) — 0, and 77 (i?) is an arbitrary magnitude, this equation 
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The second variation 5^E is of the form 



M 

5^E = J 2n {r],'n')dm, 



(57) 



(58) 



where 2rt (rj, rj') — Lr'r'Tj'^ + '2,Lrr'iiri' -\- LrrTl^ , rj' = drj/dm and lower indices denote partial derivatives with respect 
to r and r' . The value of 5^E is positive only, if the equation 



d ( dVL\ dn _ ^ 
dm \ drj' J drj ' 



has positive proper values cr^. Boundary condition is given by 



(59) 



(60) 



m=M 



For the functions £r'r', ^rr' and Lrr we find 

Lrr' — Svrrp 



du 
dp 



du 

dp/ s. 



(61) 



2p 



du 



1 f du 
P 



dg 
dr 



Using dill) and equality |Cox 1980] , [Shapiro 1983] 



9u\ [dp 



c>p^ J s^'^^ydp) s \9p)s ^ p' 



(62) 



where F is the usual generalized adiabatic exponent [Ledoux 1958j . that will be identified here with F in the 
equation of state 



{K is a constant), we obtain 



p = Kp 
Cr'r' = IGn^r'^ p'^Tp, 
Lrr' = Snrp (F — 1) , 



c„ = f r - i V — 



(63) 



(64) 



Substituting these magnitudes in (j59p . we obtain following equation for dimensionless functions ^ = rj/r: 

Q(0=^'e, (65) 



where 



m = -^^{ rpr'^) -l^m -2)p]+ pr^\ c 



1 d 

pr'^ dr \ dr J pr \ dr 



dr J 



(66) 



The operator Q is selfconjugate in the Hilbert space. Its proper values <jf (i = 0, 1 . . .) form an infinitely 
discrete sequence and are orthogonal with the weight pr'^ |Cox 1980] , Consequently, if £,k are the proper function 
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M 

1 



'yl = fJCu{Q^k)r'dm, (67) 





where ^* denotes the complex conjugate magnitude and 

M 



Jk^ J likl^r^dm. (68) 



The numerator of eq. (p7|) consist of two integral. The first of them > and the second is given by 



M 



iar(-<?) 



(3r - 4) 



4/3 



rdm. 



(69) 



Since the magnitude —g and the second term in eq. (l69|) are positive, at F > 4/3 the integral is positive. For that 
reason the proper values cr^ > 0. Thus, at least at F > 4/3, the configurations are stable. 



3.1. Numerical Results 

For a numerical solution of eq. (j65p we will rewrite it for F = Const in the dimcnsionlcss form: 



R 

X Arj 



4Ap 
R 



d^ 

dx 



n^xGM 
i?2F g 



3F — 2 x dg 
F T g dx 



(70) 



where x — r/R, Xp ~ —p/pg, ~ a{R^ /GMY/"^ is dimensionless angular frequency, and the function p{x) is the 
solution of the equation of hydrostatic equation (|55p . 

Since Xp/R ^ at a; — > 1, the magnitude in the figure brackets must tend to zero when x — > 1. Therefore, on 
the surface of the object the following equation is valid 



n^GM 3F - 2 dg/dx 



R^T g 



F 



(71) 



x = l 

Besides that, in the case of spherical symmetry, ^(0) = and we set ^(1) — 1. 

Consider, for example, the configuration with mass M = 2.6- lO^M©. Setting in eq. ([63| F = 5/3 and K = 10^^ 
(in CGS units), we find by numerical methods for eq. (j70p the three first proper functions ^(x) that arc shown in 
Fig. El 

The periods of these oscillations are given by: 



Ti = 51s, Ti = 14s, T3 = 6.7s. 



(72) 
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